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I. INTRODUCTION 



There is a growing interest in the so called analog mod- 
els, i.e. condensed matter systems which, under some cir- 
cumstances, behave as gravitational ones. The common 
feature of these systems is that the fluctuations on top of 
the ground state propagate like fields moving on a ficti- 
tious curved spacetime, whose metric tensor depends on 
the characteristics of the underlying medium. This anal- 
ogy holds both at the classical and at the quantum level. 
At the moment this setting offers the most realistic pos- 
sibility to test peculiar effects of Quantum Field Theory 
(QFT) in curved space which are still almost unobserv- 
able within the gravitational context. Among these, the 
most popular one is surely Hawking's black hole radiation 
[l|. To this aim the most promising candidates, among 
analog models, are the Bose Einstein condensates (BEC) 
@jIH> because for them the actual experimental limits are 
very close to those necessary to observe such effects Q, 
and the huge recent technological improvement in han- 
dling BEC makes us confident that this limit will be soon 
overtaken. 

In this paper we use this analogy as a new tool in an- 
alytical calculations of the quantum depletion in dilute 
weakly interacting non homogeneous BEC. The depletion 
is the number density of the non condensed particles. In 
fact, because of quantum fluctuations, even at zero tem- 
perature a non vanishing fraction of particles are non 
condensed. This quantity has been analytically calcu- 
lated for homogeneous BEC only [j| as function of the 
(constant) density of the condensate uq. The most reli- 
able attempt in extending the calculation to non homoge- 
neous BEC consists in considering the condensate locally 
as almost homogeneous, then replacing no in the result 
obtained in [j| by the actual non constant density n. This 
is the so called local density approximation (LDA) (see for 
ex. Q), which accounts for inhomogeneities only in first 
approximation, since no dependence on spatial deriva- 



tives of the density is considered. When the velocity is 
vanishing this approximation turns out to be quite good; 
however, as shown by the authors of 0, the LDA for a 
harmonic trap slightly overestimate the depletion. More- 
over, to our knowledge, nothing is known so far about ef- 
fects due to the presence of a non homogeneous velocity 
in the condensate. 

In this paper we propose a framework, borrowed from 
QFT in curved space, which allows an explicit analytical 
evaluation of how inhomogeneities and velocity affect the 
depletion in a BEC. 



II. QUANTUM DEPLETION IN BEC 

In the usual s-wave scattering approximation, a Bose 
jas is described by a field operator \I/ satisfying 
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where the coupling constant g is related to the s-wave 
scattering length a as g — 4irh 2 a/m, with m the mass of 
the constituents of the gas and V ex t the trapping external 
potential. For a dilute gas, in the mean field limit, the 
field operator ^> can be split into 



* = + tp) 



(2) 
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where VE' is the order parameter, (p the one particle fluc- 
tuation, a the distruction operator and N = ova the 
particles number operator. This splitting has been intro- 
duced following the so called number conserving approach 
HI!]. The splitting of the field # as in Eq. © leads to 
('&) = 0: this preserves the U(l) symmetry of Eq. |T]). 

As a consequence, since [0, o/V^j = 0, this framework 
keeps the total number of particles constant. However 
the use of a number conserving approach is not crucial 
for our purpouse, as will be discussed below. 
Neglecting the backreaction of fluctuations, the evolu- 
tion of the condensate is governed by the Gross-Pitaevskii 
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equation, which reads 
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where n = \* 2 \. The fluctuation part satisfies the 
Bogoliubov-de Gennes equation 



ih—ip = -— v 2 
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(4) 



The total particles density can be written as a sum of 
condensed and non condensed parts 



n to t 
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The last term describes the so called depletion. 
It is worth noticing that the depiction depends on the 
Bogoliubov-de Gennes equation only. Therefore a frame- 
work which conserves the total particles number is not 
strictly necessary (al least at first order in 1/ VN , follow- 
ing the expansion of Q). In fact the error induced by 
not using a number conserving approach is of the same 
order of that implicit in the Bogoliubov linearization of 
the theory. Instead of Eq. |(3J) one could split the field 
operator as * = * + 5*, with * = (*) and (*) = 0, 
breaking the U(l) symmetry. This will lead to the same 
results [10l | . 

For a homogeneous condensate (n = no = const) at zero 
temperature this quantum depletion can be exactly cal- 
culated [Ij]: 
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At non zero temperature one has an additional contribu- 
tion to the depletion because of thermal fluctuations. At 
low temperature (kbT <C nag, with kb the Boltzmann 
constant), summing up both contributions, one can write 
the depletion for a homogeneous condensate as 
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It is worth noting that the derivation of the zero tem- 
perature quantum depletion © is based on the full Bo- 
goliubov dispersion relation for the modes, i.e. u>k — 
ck\J\ + £ 2 fc 2 , where c = (gno/m) 1 / 2 is the velocity of 
sound and £ = ti/(^/2mc) is the healing length of the 
condensate. Unlike this, the low temperature correction 
comes from the phononic part of the spectrum only, i.e. 
LUk = ck. The integral giving the thermal depletion is 
dominated by modes with wavelength A 3> £, so the re- 
sult is insensitive to the high k behavior of the disper- 
sion relation. Therefore, while the correct description 
of the ground state requires the knowledge of the full 
"microscopic" theory deviations from the ground 

state value can, in some circumstances, be evaluated with 
sufficient accuracy using the low energy "macroscopic" 



phonon theory. This feature of the low energy theory is 
familiar in the Casimir effect. The standard Casimir en- 
ergy density between two conducting plates can be eval- 
uated with sufficient accuracy using the photon disper- 
sion relation, provided the plate separation is much larger 
than the interatomic distance of the constituents of the 
plates [13 . 

Moving to non homogeneous condensates, when the cur- 
rent j = ^^V^ — c.cj vanishes, the depletion is 

usually approximated by Eq. with now no replaced 
by the actual space-time varying density n{x, t) in the 
local density approximation (LDA) [f|. We indicate this 



as 
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where the first term is the zero temperature quantum 
depletion and the second the thermal depletion, both in 
LDA. 



III. BEC AND THE GRAVITATIONAL 
ANALOGY 

The aim of this paper is to show how, in view of the 
analogy described at the beginning, through the use of 
QFT in curved space and under some hypothesis, it is 
possible to obtain an estimate of the depletion which ex- 
plicitly accounts for the inhomogeneities of the conden- 
sate (i.e. din ^ 0) and for a non trivial velocity (i.e. a 
velocity that can not be simply eliminated by a Galileo 
transformation) . 

For BEC the analogy traces back to the hydrodynamical 
approximation of the basic equations for the condensate 
([31 0]) . In the density-phase representation for the con- 
densate field \P = ^/ne l8 / R , Eq. © is split in a continuity 



and an Euler equation, namely 
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where v — VO/m is the velocity field. The last term in 
(fl~0|) is called quantum pressure. The hydrodynamical 
approximation consists in neglecting this compared to 
the mean field term gn. For this approximation to be 
valid, the typical length scale L for the variations of the 
condensate density n has to be much bigger than the 
healing length (L ^> £). Decomposing the fluctuation 
operator ip as 



<p = e ie ' h 



2y/n n 
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where hi and 9\ are real fields, one finds that for wave- 
lengths A 3> £, the Bogoliubov-de Gennes Eq. Q can be 
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written as 



□0i = , 
h x = -g- 1 (dA + m- x V0 • V0i 



(12) 
(13) 

where □ = - 7 =d M (-\/ = p(7' u/ 9 ! /) is the covariant 
d'Alembertian calculated from the so called "acoustic 
metric" 



SV = — 

mc 



— (c 2 — V ) —v 
-v 1 



(14) 



where c is the local sound speed related to the density by 
mc 2 = gn. Hence the fluctuation 0i behaves exactly as 
a massless minimally coupled scalar held propagating on 
a fictitious curved spacetime described by the acoustic 
metric g^, whose line element reads 

ds 2 = — [-c 2 dt 2 + (dx - vdt) ■ (dx - vdt)] . (15) 
mc 

This is the core of the analogy [HI]. 
For a homogeneous condensate with constant velocity vq, 
the acoustic spacetime is flat as can be seen performing 
the Galileo transformation dx — ► dx — vodt. This conden- 
sate configuration is a sort of ground state of the theory. 
On the other hand inhomogeneities or a non trivial ve- 
locity produce curvature. 

We will focus our analysis on the latter kind of configu- 
rations, looking for the corrections on the depletion in- 
duced by curvature with respect to the homogeneous (let 
say flat) result. 



IV. QFT IN CURVED SPACE AND THE 
DEPLETION IN BEC 



As just mentioned, phonons in a BEC propagate as a 
massless scalar field on a curved background spacetime: 
therefore they can be described through the gravitational 
formalism. In the case of non homogeneous BEC, the 
curvature is not zero and we expect modifications with 
respect to the standard result induced by this curvature. 

In this case modes with A < £ <C L see the condensate 
as almost homogeneous. After renormalization, they will 
give the quantum depletion term in LDA, n q LDA . De- 
viations from this come from modes with wavelength 

A > L ^s> £. These are indeed the modes who feel the 
inhomogeneities of the condensate. They behave as a 
massles scalar field on the curved spacetime g^ v . So the 
low energy phonon theory (even if unable to reproduce 
the ground state value of the depletion) can describe with 
sufficient accuracy how inhomogeneities influence the de- 
pletion. Now, within this spirit, calling h v h the depletion 
induced by the phonons, we can write h = n q LDA + n p h, 
where, using Eq. (fTTj) : 



n-i/ph 
An 
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and the expectation values are calculated using the 
phonon theory. The last approximation in Eq. (1161) 
comes from taking h\ from Eq. (|13p and estimating 

V0i « | and dA « -if- c - 

In this way the calculation of the depletion has been re- 
duced to the calculation of the expectation value of the 
squared of a massless scalar field minimally coupled to 
the acoustic spacetime (IT4")) . n(9 2 ) p h/h 2 is dominated by 
wavelengths ~ L, therefore one can safely extrapolate the 
linear dispersion relation to high frequencies without af- 
fecting the finite rcnormalized result, since the diverging 
terms appearing are then subtracted away. 
In QFT in curved space a useful method of renormaliza- 
tion is point-splitting [l4j. Within this procedure, (8f(x)) 
(here and in the following we will drop the subscript ph 
from the expectation values) is the coincidence limit of 
the Euclidean Green function 



Re ( lim Ge(x, x') 

\x f — *x 



(17) 



Renormalization is performed by subtraction of the diver- 
gent part, the so called De Witt-Schwinger term Gds, 
obtained by the short-distance expansion of the Green 
function: 

0l(x)) ren = lim {{el{x,x')) unren - G DS {x,x')) . (18) 

X — *x 

Gds{x,x') is purely geometrical, state independent and 
contains all the short-distance divergences of the theory 
(for technical details see for ex. [Hj]). 



A. Static spherically symmetric configuration 

We now restrict our discussion to a static spherically 
symmetric configuration (n = n(r), v = 0). The resulting 
acoustic metric is invariant under time-translation and 
dt is the Killing vector associated to this symmetry. The 
acoustic metric (|14[) can be put in the standard form 



ds 2 = -f{R)dt 2 + h(R)dR 2 + R 2 dn 2 , (19) 
by the coordinate transformation 

(20) 



R=J— r. 
mc 



Note that this coordinate transformation is valid only if 
R 7^ const. The R = const case will be treated sepa- 
rately. Now 



m 2c 



(21) 



r, 2 



and ' = d r . 

For configurations of this kind, it is natural to consider 
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thermal equilibrium states, i.e. whose Euclidean Green 
function is periodic in T _1 in the imaginary time. For 
these states, Anderson et al. [HI have developed a WKB 
technique to evaluate (Ol(x,x')) unren in the case of a 
generic metric of type (|19l) leading to 



WKB 

ren 



,T 2 



12ft/ 



11 
96^r 



/'ft' 



Hv 2 f) 



f" 



f 



'2 



96Tr 2 f 2 h 



r 



192n 2 fh 2 96n 2 fh A%it 2 Rfh 



(22) 



where only here ' = 5r, while 1Z is the curvature scalar 
and n an arbitrary scale introduced by renormalization 
[l4| . The first term describes a gas of massless scalar par- 
ticles in thermal equilibrium at the temperature T in a 
curved spacetime. The factor 1// accounts for the "grav- 
itational red-shift" (Tolman law [IB]). The other terms 
describe how the quantum fluctuations are affected by 
the curvature of the effective spacetime. Entering our ex- 
pressions (|21[) in Eq. ([22]) one gets the following analytic 
expression for the thermal and inhomogeneity dependent 
part of the depletion: 



n p h = 
n 2 B T 2 



m 
gn 



12ft 3 
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+ 487T 2 ft 
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l/2 n 3/2 
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Note that the first term in Eq. 



elegantly reproduces 



l lda- 



The other terms give the corrections to the depletion in 
LDA induced explicitely by the spatial derivatives of the 
condensate density. They all vanish for a homogeneous 
BEC. 

One can obtain a numerical estimate of the relative 
weight of the inhomoheneities induced terms. For N = 
10 6 Rb atoms trapped by a harmonic potential (L osc = 
lum, a = 5nm, Ro = L{\haN / L) 1 /^ the size of the con- 
densate) Eq. (|23|) gives at T = a correction with respect 
to the LDA result of < 0.1% at the centre of the trap, 
reaching ~ 1% approaching the boundary (r ~ 0.9i?o)- 
For N — 10 3 atoms the correction will be almost 2% 
at the center of the trap and almost 20% at r ~ 0.7i?o 
[H| ■ In both cases the corrections turn out to be neg- 
ative, therefore lowering the actual number of non con- 
densed particles. This is in qualitative agreement with 
what pointed out in 



B. Radial velocity 



In the presence of a nonvanishing radial velocity (i.e. 
= u(r)X) pjl the system is no longer in equilibrium 



and no analytical expression for {9 2 ) ren is known. One 
can only evaluate how vacuum fluctuations are modified 
by the presence of the radial velocity. Being the den- 
sity and the velocity functions of the radial coordinate 
r only, one can transform the acoustic metric (|14p in 
the standard spherically symmetric form of Eq. (fill)) by 
the following coordinate transformation (assuming again 
R 7^ const) 



T = t + 



v dr 



R=\—r, (24) 



leading to a diagonal form of the type of Eq. (|19p with 
now 

„2 J 



I 



n 



mc 
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(25) 



The vacuum part of Eq. (|22[) for this new configuration 
gives 

"3 
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p' 1 96ir 2 hc 2 
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(3c 2 c' ~ 2v 2 c' - cvv') - —cc' 2 ~ -v 2 c r - 
4 4c 

+ l c 2 c" - \ 2 C " - cv ' 2 - cvv "+ 
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—vc'v' + 



+ i(3c 2 c'-f 2 c'-2cw')| . (26) 

This equation generalizes n p h of Eq. (|23p at T = in 
the presence of a radial velocity. Note that for homo- 
geneous condensate (c = const) this expression does not 
vanish for |tj| = const, as a consequence of the fact that a 
radial velocity can not be eliminated by a Galileo trans- 
formation. Eq. |J2g} reduces to ||23J) (at T = 0) only for 
v(r) = 0. 

It is important to stress that Eq. (|26[) is strictly valid 
only in regions where the fluid velocity does not reach 
the sound velocity. For |u| = c, Eq. (f2"6")l diverges. This 
local divergence has nothing to do with the ultraviolet di- 
vergences already canceled out in deriving Eq. (|26p . This 
divergence is quite familiar in the gravitational context: 
in fact for \v\ = c the goo term in the acoustic metric 
pT|) vanishes and this is a hint of the fact that a hori- 
zon forms (the so called sonic horizon). This reflects on 
the modes of the scalar field that turn out to be singular 
on the horizon 20]. As a direct consequence, a singular 
vacuum polarization on the horizon of a black hole ap- 
pears. However it is well known that taking the black 
hole formation process into account a Hawking radiation 
emerges: its presence is responsible for an additional con- 
tribution to (9 2 ) U nren which has an equal and opposite 
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sign divergence, leading to a (9\) ren regular on the hori- 
zon. This has been explicitely shown for a Schwarzschild 
black hole (for which f(r) = g' 1 ^) = 1 - 2M/r) [H]. A 
similar behavior is expected also in the BEC case. 
Unfortunately for a general metric as in Eq. (fT9|) an ex- 
plicit analytic form for the {0 2 } unren taking Hawking ra- 
diation into account and valid in the near horizon region 
is still unknown. 
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For the sake of completeness let us consider the case 
when the coordinate transformation (|20p is singular, i.e. 
when n/c oc 1/r 2 . The acoustic metric (fT9"|) shows that 
the surface (r, t) — const has constant area. In this case 
Eq. ([22]) is replaced by 
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967T 2 /ft_ 

2 ) _1 and / unchanged as in 



with ' = d r , ft : 
Eq. (EH]). 

For the standard equation of state for BEC mc 2 = gn. 
one can take n = n^r^/r^, with no and ro two constants 
having the dimension of density and length respectively. 
The resulting quantum depletion reads 



48ir 2 hr 4 



(n gm) 
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V. CONCLUSIONS 

Summarizing, we have shown how sophisticated renor- 
malization techniques originally developed to study 
quantum effects in black hole physics, can be ex- 
tremely useful in a (seemingly) complete different con- 
text, namely BEC. In this paper these techniques have 
been used to give an analytical estimation of the effects of 
spatial derivatives and a non trivial velocity on the deple- 
tion of a BEC. Our predictions are beyond the present ex- 
perimental capability. Nonetheless this work represents 
the first attempt to go beyond the LDA, providing a gen- 
eral analytical scheme that can be used in cold-atoms 
systems within the hydrodynamical approximation. 
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